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SUMMARY 


The  interaction  of  electromagnetic  fields  with  biological  tissues  is  of  increasing 
importance  from  the  standpoint  of  potential  health  hazards  and  of  possible  beneficial 
(e.g.,  diagnostic  and  therapeutic)  effects.  To  ".aoress  these  issues,  the  research  ac¬ 
complished  during  this  period  focused  on  the  biosynthesis,  assembly,  and  catabolism 
of  highly  specialized  protein  and  polysaccharide  molecules  by  mammalian  connec¬ 
tive  tissue  cells  the  synthesis  and  processing  of  these  macromolecules  is  one 

of  the  most  important  biological  functions  of  connective  tissue  cells,  altered  matrix 
metabolism  is  known  to  be  a  most  sensitive  indicator  of  the  effectiveness  of  an  envi¬ 
ronmental  stimulus,  and  of  the  pathway  by  which  the  stimulus  evokes  a  response. 

The  long-term  goal  is  to  study  long-term  (chronic)  and  short-term  (acute)  ef¬ 
fects  of  electromagnetic  field  exposure  on  the  rynthesis  and  turnover  of  highly  charged 
proteoglycan  (PG)  molecules  and  their  glycosaminoglycan  (GAG)  constituents,  col¬ 
lagens,  and  other  non-collagenous  proteins  by  normal  cartilage  cells  (chondrocytes) 
extracted  from  articular  cartilage,  and  in  rat  chondrosarcoma  cells  (a  continuous  cell 
line).  Another  long-term  goal  is  to  study  the  possible  role  of  fields  in  stimulating 
growth  factor  synthesis,  and  the  combined  role  of  fields  and  growrth  factors  on  chon¬ 
drocyte  metabolism. 

To  address  these  goals,  the  specific  objectives  of  this  research  period  were:  (1) 
determine  whether  chondrocytes  can  be  maintained  in  long-term  agarose  culture,  in 
order  to  allow  experiments  involving  chronic  exposure  to  fields;  (2)  quantify  chon¬ 
drocyte  proliferation,  proteoglycan  and  protein  biosynthetic  rates,  total  proteoglycan 
accumulation,  and  matrix  physical  properties  during  long-term  culture;  (3)  construct 
a  special  chamber  to  enable  application  of  electric  fields  to  agarose/chondrocyte  disks 
maintained  in  culture  in  a  standard  incubator;  and  (4)  begin  electrical  stimulation 
experiments  over  a  range  of  current  density  amplitudes. 


INTERACTION  OF  ELECTROMAGNETIC  FIELDS  WITH 
BIOSYNTHETIC  PROCESSES  IN  CONNECTIVE  TISSUE  CELLS 


INTRODUCTION 

Articular  Cartilage 

Articular  cartilage  is  the  dense  connective  tissue  that  functions  as  a  bear¬ 
ing  material  in  synovial  joints.  Adult  articular  cartilage  is  avascular,  aneural,  and 
alymphatic;  cell  nutrition  is  derived  primarily  from  the  synovial  fluid  [1].  The  chon¬ 
drocytes  are  responsible  for  the  synthesis,  maintenance,  and  gradual  turnover  of  an 
extracellular  matrix  (ECM)  composed  principally  of  a  hydrated  collagen  fibril  net¬ 
work  enmeshed  in  a  gel  of  highly  charged  proteoglycan  molecules  (PG).  The  compo¬ 
sition  and  architecture  of  the  matrix  [2]  and  the  tissue’s  high  water  content  (70-80% 
of  wet  weight)  enable  cartilage  to  withstand  complex  compressive,  tensile  and  shear 
deformation  in  joints  [3,  4]. 

Cartilage  ECM  can  remodel  to  meet  the  functional  demands  of  loading.  It 
is  known,  for  example,  that  electrostatic  and  osmotic  interactions  between  the  pro¬ 
teoglycan  constituents  of  the  ECM  enable  cartilage  to  resist  compressive  loads  [2,  5]. 
Animal  studies  have  shown  that  PG  content  is  higher  in  cartilage  that  is  habitually 
loaded  [6],  while  immobilization  of  joints  leads  to  a  decrease  in  PG  synthesis  and  total 
PG  content  [7],  Osteoarthritic  cartilage  initially  shows  enhanced  matrix  synthesis  in 
an  apparent  attempt  to  replace  degraded  PG  molecules  that  have  diffused  out  of  the 
tissue.  The  mechanisms  that  mediate  the  biological  response  to  mechanical  loads  are 
not  understood. 

Biosynthesis,  metabolism,  and  turnover  of  cartilage  matrix  components  have 
been  studied  extensively  in  vivo  and  using  cell  and  tissue  culture  models  [8].  Recent  in 
vitro  studies  have  shown  that  mechanical  loads  can  also  induce  metabolic  responses  in 
cartilage  organ  cultures.  Static  compression  has  been  found  to  decrease  proteoglycan 
and  protein  synthesis  [9,  19,  11],  while  dynamic  compression  at  certain  frequencies 
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and  amplitudes  stimulated  synthesis  of  these  matrix  constituents  [11].  Taken  together, 
these  data  suggest  that  the  presence  of  physiologic  mechanical  loading  forces  may  be 
necessary  for  the  proper  assembly  and  maintenance  of  a  matrix  that  can  function  in 
the  joint. 

Chondrocytes  Cultured  in  Agarose 

Investigators  have  recently  found  that  bovine  [12,  13,  14],  rabbit  [15],  and 
human  [16]  chondrocytes  as  well  as  Swarm  rat  chondrosarcoma  cells  [17]  isolated 
aid  cultured  in  agarose  gel  maintain  their  phenotype  as  evidenced  by  the  synthe¬ 
sis  of  type  II  collagen  and  large  aggregating  proteoglycans.  The  newly  synthesized 
macromolecules  are  effectively  trapped  within  the  gel.  When  placed  in  monolayer  cul¬ 
ture,  rabbit  and  human  chondrocytes,  and  Swarm  rat  chondrosarcoma  cells  attach, 
flatten,  spread,  and  dedifrerentiate  by  beginning  to  synthesize  type  I  collagen  and 
small  non-aggregating  proteoglycans.  Matrix  macromolecules  produced  in  monolayer 
culture  diffuse  out  into  the  medium  and  are  not  retained  to  form  an  intact  matrix. 
(Although  cell  phenotype  may  be  maintained  in  monolayer  culture  for  bovine  chon¬ 
drocytes,  a  rounded  cell  shape  is  not  [18].)  Cell  shape  has  been  implicated  as  the 
primary  determinant  of  the  phenotype  and  recent  studies  [19,  2.0]  have  suggested 
that  the  integrity  of  the  microfilament  cytoskcletal  structure  may  be  more  proximal 
to  phenotypic  change  than  is  cell  shape. 

Agarose  culture  of  chondrocytes  has  been  used  as  an  in  vitro  model  for  studying 
matrix  degradation  (induced  by  retinol)  which  may  resemble  that  which  occurs  in  vivo 
in  early  stages  of  osteoarthritis.  Possible  therapeutic  agents  such  as  anti-invasion 
factor,  Arteparon,  Eglin  C,  and  pepstatin  (most  are  primarily  protease  inhibitors) 
have  been  tested  [12].  The  ability  to  selectively  seed  certain  cell  populations  has 
allowed  the  study  of  differences  between  cell  lubpopulations.  Cells  from  the  superficial 
zone  (top  20-40  fim)  of  bovine  metacarpal  ohalangea!  cartilage  were  shown  to  have  a 
more  irregular  shape,  to  divide  slower,  to  synthesize  and  retain  less  proteoglycan,  to 
synthesize  a  smaller  and  less  aggregable  proteoglycan,  and  to  produce  larger  collagen 
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fibrils,  than  chondrocytes  from  deeper  zones  [13,  14).  Transparency  of  the  agarose 
system  has  allowed  the  study  of  cell  size  over  the  culture  period  [17]  and  the  study 
of  colony  formation  in  response  to  insulin-like  growth  factor-I  and  growth  hormone 
[21].  Monolayer  culture  has  been  used  to  amplify  cell  number  prior  to  placement  in 
agarose  and  subsequent  phenotype  reexpression  [16]. 

To  date,  most  studies  have  used  cell  densities  10-100  times  less  than  that  found 
in  native  cartilage.  Although  the  synthesis  and  accumulation  of  normal  matrix  com¬ 
ponents  has  been  observed  biochemically,  no  attempt  has  been  made  to  measure  the 
mechanical,  electromechanical,  and  physicochemical  properties  of  the  matrix.  Nev¬ 
ertheless,  these  studies  suggest  that  the  agarose/chondrocyte  system  may  be  ideally 
suited  for  examining  the  synthesis  of  a  mechanically  functional  cartilage-like  matrix. 

The  agarose/ chondrocyte  system  may  also  benefit  the  examination  of  trans¬ 
duction  mechanisms  which  mediate  cell  response  to  mechan. cal  and  electrical  stimuli 
to  the  tissue.  The  amount  of  matrix  present  during  stimulation  may  be  controlled 
by  choosing  different  time  points  during  culture.  For  example,  the  choice  of  an  early 
time  point  will  highlight  cell  compression  and  reduce  physicochemical  and  direct  cell- 
matrix  interactions.  Transparency  of  the  system  may  allow  measurements  of  single 
cell  properties  such  as  cell  shape  or  those  exhibited  by  fluorescent  probes.  And  lastly, 
the  capacity  to  control  the  population  of  cells  in  a  given  experiment  will  have  numer¬ 
ous  applications. 

Continuum  Models  of  the  Physical  Properties  of  Cartilage  Matrix 

The  physical  properties  of  cartilage,  that  is  its  large  compressive  stiffness, 
the  existence  of  electrokinetic  phenomena  and  electrostatically  influenced  internal 
swelling  pressures  and  ion  distributions,  can  be  thought  of  as  resulting  from  a  highly 
charged  proteoglycan  gel  enmeshed  within  a  collagen  framework.  The  PG,  containing 
a  large  number  of  fixed  charge  groups  which  repel  each  other,  has  a  tendency  to  swell 
in  the  native  state  but  is  contained  by  tensile  stress  within  the  collagen  fibril  network. 
Although  the  collagen  is  vital  in  not  only  containing  the  PG’s  but  in  resisting  tensile 
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and  shear  strain,  it  is  the  PG  that  is  the  determinant  of  many  properties  related 
to  compressive  stress.  In  the  physiologic  state,  electrostatic  repulsion  between  fixed 
charge  groups  contributes  about  one-half  of  the  equilibrium  modulus  [22].  Transport 
properties  such  as  hydraulic  permeability,  electrical  conductivity,  and  electrokinetic 
coefficients,  are  also  significantly  affected  by  interactions  within  the  proteoglycan 
component  of  the  matrix  [23,  24], 

Macrocontinuum  Theories 

Most  theories  pertaining  to  the  mechanical,  electromechanical,  and  physico¬ 
chemical  properties  of  cartilage  matrix,  and  other  charged,  porous  materials,  have 
been  macrocontinuum  theories.  Such  theories  inherently  ignore  molecular  structure 
by  averaging  over  a  volume  element  large  enough  to  contain  many  molecular  com¬ 
ponents.  The  characteristic  dimension  is  then  much  larger  than  the  Debye  length 
(~  8A°  in  physiological  saline)  so  that  quasineutrality  is  assumed  valid  everywhere 
except  over  a  small  region  at  the  tissue- liquid  interface  [5]. 

Equilibrium  properties,  such  as  ionic  distributions  and  osmotic  swelling  pres¬ 
sure,  are  characterized  by  Donnan  equilibrium  in  which  a  spatially  invariant  electro¬ 
static  potential  exists  (due  to  fixed  charge)  in  the  tissue  and  is  different  from  that 
in  the  bath  [25].  The  exponential  nature  of  the  Boitzmann  distribution  demands 
that  more  counterions  are  attracted  than  coins  are  repelled  resulting  in  a  higher  total 
ionic  concentration  in  the  tissue  than  in  the  bath.  The  excess  ionic  concentration  is 
interpreted  as  a  Donnan  osmotic  swelling  pressure  using  van’t  Hoff’s  relation.  Mathe¬ 
matically  the  model  is  very  simple,  being  two  algebraic  relations  requiring  knowledge 
of  bath  ionic  strength,  the  fixed  charge  density,  and  mean  ionic  activity  coefficients 
both  in  the  tissue  and  in  the  bath.  The  disadvantage  is  that  equilibrium  measurements 
of  the  activity  coefficient  in  the  tissue  requires  the  assumption  of  Donnan  theory  so 
that  the  validity  of  the  theory  itself  is  not  determined  [26,  27,  28,  29].  Some  attempts 
have  been  made  to  derive  expressions  for  the  activity  coefficients  based  on  molecular 
structure  using  polyelectrolyte  condensation  theory  [30,  31].  These  attempts  have 
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met  with  limited  success  and  are  excursions  into  a  microcontinuum  description 
Microcontinuum  Theories 

The  effort  in  the  microcontinuum  approach  is  to  relate  properties  to  molecular 
structure.  Judicious  choice  of  the  molecular  model  is  important  in  order  to  obtain  a 
solvable  set  of  equations  that  still  captures  the  elements  of  molecular  structure  that 
give  rise  to  particular  properties.  Given  the  predominance  of  the  proteoglycan  as 
the  main  contributor  to  osmotic  swelling  pressure,  ion  distribution,  and  transport 
properties,  the  choice  to  be  made  is  one  of  geometry. 

The  basic  molecular  unit  of  the  PG  is  the  chondroitin  sulfate  (a  glycosamino- 
glycan  or  GAG)  disacharide  chain.  Each  disacharide  contains  two  charges,  a  carboxyl 
and  a  sulfate,  both  ionized  at  physiological  pH.  Perhaps  50  chains  of  about  40  dis- 
acharides  each  are  covalently  linked  to  a  core  protein  forming  a  brush-like  proteogly¬ 
can  monomer  of  —  108  MW  which  is  secreted  by  the  cell  as  a  unit.  To  contribute 
functionally  to  tissue  properties,  the  monomer  must  find  a  hyaluronic  acid  chain  to 
which  a  noncovalent  bond  is  then  formed,  being  stabilized  by  a  link  protein,  to  form 
a  ^ro'^oglycan  aggregate  i32l.  The  proteoglycan  is  then  effectively  trapped  within 
the  collagen  network  and  can  not  diffuse  out  of  the  tissue.  A  geometrical  model  of 
an  aggregate,  a  monomer  or  even  a  finite  length  chondroitin  sulfate  chain  would  pro¬ 
duce  considerable  mathematical  difficulties  in  solving  the  relevant  Poisson- Boltzmann 
equation.  The  highest  level  of  tractable  complexity  is  that  of  a  rod-like  polyelectrolyte, 
representing  a  cross  section  of  a  chondroitin  sulfate  ( C S )  chain. 

The  CS  chain  is  of  the  correct  length  scale  when  considering  properties  arising 
from  electrostatic  repulsion.  Electrostatic  forces  between  charged  bodies  in  an  isotonic 
electrolyte  solution  in  equilibrium  are  completely  shielded  across  distances  greater 
than  a  few  nanometers  (CS  diameter  ~  1  nm).  If  the  rod  is  considered  to  be  a 
monomer,  with  a  diameter  of  ~  100  nm,  any  reasonable  degree  of  hydration  would 
place  monomers  far  enough  apart  to  shield  any  electrostatic  interactions.  This  is  not 
the  case  in  the  tissue.  Electron  microscopy  of  intact  tissue  has  indicated  that  PG 


monomers  are  packed  so  tightly  that  a  CS  chain  on  one  monomer  may  be  as  close  to 
a  CS  chain  on  another  monomer  as  it  is  to  an  adjacent  chain  on  the  same  monomer 
[33]. 

The  microscopic  picture  is  then  a  meshwork  of  CS  chains  stabilized  in  place  by 
the  hierarchy  of  monomer  (core  protein),  aggregate  (hyaluronic  acid),  and  collagen 
interactions.  There  may  be  some  non-random  orientation  of  the  meshwork  giving  rise 
to  periodic  arrays  or  nematic  phases  [34],  at  least  within  limited  domains.  However, 
as  a  first  approach,  the  assumption  of  an  isotropic  phase  is  accompanied  by  sufficient 
complexity.  The  model  is  then  the  cell-model  where  a  charged  rod  of  infinite  extent 
is  surrounded  by  an  annulus  of  electrolyte.  The  ratio  of  the  two  radii  is  determined 
by  polyelectrolyte  concentration,  or  equivalently  by  tissue  hydration. 

CHONDROCYTES  IN  AGAROSE 

Summary 

We  have  studied  the  ability  of  cartilage  cells  extracted  from  native  cartilage 
to  synthesize  and  accumulate  a  normal,  cartilage-like  extracellular  matrix  in  agarose 
gel  culture.  We  have  characterized  the  extent  and  the  time  evolution  of  chondro¬ 
cyte  proliferation,  synthesis  of  GAG  and  proteins,  loss  of  GAG,  and  total  deposition 
of  G AG-containing  matrix  within  agarose  gels  during  long-term  culture.  To  assess 
whether  the  matrix  deposited  within  the  agarose  gel  is  a  mechanically  functional 
cartilage-like  matrix,  we  measured  several  important  mechanical  and  electromechan¬ 
ical  properties  of  the  agarose/ chondrocyte  disks  at  selected  times  during  long-term 
culture:  equilibrium  elastic  modulus,  dynamic  stiffness,  and  streaming  potential  in¬ 
duced  by  oscillatory  mechanical  compression.  The  results  of  these  studies  suggest 
that  (1)  both  normal  chondrocytes  and  Swarm  rat  chondrosarcoma  cells  in  agarose 
culture  can,  under  proper  culture  conditions,  continue  to  synthesize  matrix  macro- 
molecules  at  a  rate  similar  to  that  in  native  cartilage,  and  (2)  chondrocytes  in  agarose 
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can  successfully  mediate  the  assembly  and  accumulation  of  a  normal,  mechanically 
functional  extracellular  matrix. 

Previous  experiments  of  chondrocytes  or  chondrosarcoma  cells  cultured  in 
agarose  utilizing  a  cell  density  ranging  from  (1  -  20)  X  10flcellc/ml  and  an  agarose 
concentration  in  the  range  of  (1  —  5)%  cultured  in  serum  supplemented  Dulbecco’s 
Modified  Eagles  Medium  (DMEM)  have  identified  an  optimal  system  consisting  of 
20  x  l06cclls/ml  in  2%  agarose.  Two  experiments  using  normal  chondrocytes  with 
these  parameters  have  been  completed  yielding  similar  results.  Experiments  using 
Swarm  rat  chondrosarcoma  cells  have  also  been  completed. 

Methods 

Articular  cartilage  from  the  femoropatellar  groove  from  1-to  2-week-old  calves 
was  harvested  as  previously  described  [35]  and  incubated  in  DMEM  supplemented 
with  12.5  mM  Hepes,  0.1  mM  nonessential  amino  acids,  0.4  raM  L-proline,  10% 
FBS,  50  /i g/ml  ascorbate,  and  0.1%  penicillin /streptomycin  changed  daily  at  5% 
C02,  37° C.  One  week  later  cells  were  extracted  by  sequential  pronase  and  collagenase 
digestion  [18].  In  addition,  Swr  ,m  rat  chondrocytes  were  prepared  for  incorporation 
into  agarose  gels  as  previously  described  [17]. 

Chondrocytes  were  mixed  into  media  containing  2%  (w:v)  low  melting  tem¬ 
perature  (FMC)  SeaPlaque  agarose  [12,  17]  at  ~  2  x  107cells/ml  and  cast  at  37°C 
between  slab  gel  electrophoresis  plates  separated  by  1-mm-thick  Teflon  spacers,  as 
shown  schematically  in  Figure  1.  After  gelling  at  4°C  for  2  h,  16-mm-diameter  by 
1-mm-thick  disks  were  cored  from  the  gel.  The  chondrocyte/agarose  disks  were  sub¬ 
sequently  cultured  on  top  of  nylon  mesh  (to  promote  nutrient  diffusion  from  below). 
Media  (as  above)  was  changed  daily  and  analyzed  for  GAG  content  by  dimethyl- 
methylene  blue  (DMB)  dye  binding  [36],  Control  disks  without  chondrocytes  were 
prepared  and  maintained  in  the  same  manner. 

The  mechanical  and  electromechanical  properties  of  chondrocyte/ agarose  disks 
and  control  disks  were  measured  in  uniaxial  confined  compression  using  the  test  cham- 
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Preparation  of  Chondrocyte/Agarose  Plugs 

~2  Days 


Calf  Articular  Cartilage 

femoropatellar  groove 


Cells  in  Gels  Casting  Apparatus 

(cross-section) 


Figure  1.  Preparation  of  chondrocyte/agarose  cultures. 
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Mechanical  and 
Electromechanical 
Testing 


Biochemical  Analysis 


12  x  3mm  plugs 


•  •  •  • 

•  •  •  •  ^ 

•  •  •  • 


•DNA  content 
•GAG  content 
•Proline  incorp. 
•Sulfate  incorp. 
•Alcian  Blue 


Figure  2.  Mechanical,  electromechanical  and  biochemical  tests. 


-  n  senematiraliy  in  Figure  2,  which  was  mounted  in  a  Dynastat  mechanical 
-  t  Several  i3-mm  diameter  disks  were  cored  from  the  larger  disks  re- 
...  '  j i-r  at  selected  intervals,  and  placed  between  silver/silver  chloride 

*  .••.des  separated  by  a  porous  platen,  all  bathed  in  phosphate-buffered 
’  -  f  ’  A  20  A  static  offset  strain  was  applied  and  the  load  recorded 

•  >.-  rel  '...x  at  ion  reached  equilibrium.  This  equilibrium  load  and  strain  was  used 
-‘mipute  the  equilibrium  elastic  confined  compression  modulus.  A  0.5%  sinusoidal 
strain  was  then  superimposed  at  frequencies  between  0.01  and  1  Hz,  and  the  result¬ 
ing  oscillatory  load  and  streaming  potential  were  recorded.  Disks  were  then  removed, 
weighed,  lyophilized,  and  frozen  and  reweighed.  The  hydraulic  permeability,  elec- 
f rckinetic  coefficient,  oscillatory  stiffness,  and  streaming  potential  were  computed 
f:.,m  the  data  as  described  previously  [37,  23], 

Each  week,  groups  of  9  to  12  3-mm-diameter  disks  were  separately  incubated 
for  16  h  in  media  containing  10  /xCi/ml  [35Slsulfate  and  20  /iCi/ml  [5-3H]proline  to 
assess  GAG  and  protein  synthesis,  respectively  (Fig.  2).  These  disks  were  washed  in 
l)A;rM  and  digested  with  papain.  Radiolabel  incorporation  was  assessed  by  scintil- 
r.m  OKing,  and  GAG  was  assessed  by  the  DMB  dye  binding  assay.  In  addition, 
m  rit  of  deoxyribonucleic  acid  (DNA)  was  measured  in  these  same  specimens 
!  :  sf  a  fluorescence  enhancement  assay  using  the  Hoechst  33258  dye  [38]. 

Results  and  Discussion 

.X* aural  Chondrocytes 

Table  1  shows  the  biochemical  properties  of  agarose/chondrocyte  disks  on  day 
50  of  a  long  term  culture  that  was  maintained  for  70  days.  These  properties  are 
compared  to  those  of  intact  calf  articular  cartilage  similar  in  age  and  location  to 
that  from  which  the  cells  were  extracted.  In  agarose/chondrocyte  disks,  total  GAG 
content  increased  7-fold  over  10  weeks  in  agarose  culture  (Fig.  3)  reaching  ~25%  of 
th:  1  in  native  parent  cartilage  disks  (Table  1).  By  day  30,  the  amount  of  GAG  alone 
was  approximately  the  same  as  the  agarose  content  of  agarose/chondrocyte  disks. 
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Fable  ’.  CHONDROCYTES  IN  AGAROSE  :  BIOCHEMICAL  PROPERTIES  ON 

DAY  50 


COMPARISON  : 


CHONDROCYTE/AGAROSE  DISKS 
and  CALF  CARTILAGE  EXPLANTS 


Chondrocyte/ 
Agarose  diskst 

Calf  Cartilage 
‘parent  disks’ 

Total  DNA  inc-pi  disk  vol.) 

0.14  ±  0.02  (12) 

0  56  x  0.13  (30;  [1]  ' 

Ceil  Density  ;io”celi/ml)  (7.7pg/celi(7j) 

.  20.2  ±  2.8  (12; 

73  x  17  (30)  [ij 

Total  GAG  (pg/pi  disk  voi.) 

14.1  x  2.6  (12) 

57-11  (30)  [1i 

GAG  release  rate  (pg/(pl  disk  vol)/day) 

|  0.20  ±  0.03  (12) 

0.62  x  0.12  (12)  [2( 

GAG  accumulation  rate  (pg/(pi  disk  vol)/day)  j  0.22  x  0.03  (121 

1.65  [1] 

Sulf-incorp.  (nmol  sulfate/ io6cell/day) 

22  ±4  (12) 

45  x  21  (2C)  (1] 

Prol-incorp.  mmol  proiine/iCFcell/day) 

17x2(12) 

20.1  ±  9.5  (28)  [1] 

t  data  for  day  50  MEAN  ±  SD  In) 

1  j  Sah  RL  et  al.  J  Orthop  Res,  7:619,  1989. 

21  Sah  RL  et  al.  Trans  ORS,  14:50.  1989. 


DNA  content  appeared  to  increase  by  a  factor  of  2  over  the  10-week  period  (Fig.  3). 
The  .significant  increase  in  disk  wet  weight  and  thickness  (Fig.  4)  with  time  in  culture 
is  a  direct  consequence  of  the  accumulation  of  matrix  and  the  high  osmotic  swelling 
pressure  of  the  GAG  constituents. 

Sulfate  and  proline  incorporation  (Fig.  5)  reflects  the  rate  of  synthesis  of  pro¬ 
teoglycans  and  total  protein,  respectively.  Sulfate  and  proline  incorporation  levels 
in  agarose  culture  (Fig.  5)  were  similar  to  that  in  the  intact  parent  cartilage  (Ta- 
ble  1  .  The  decrease  in  radiolabel  incorporation  after  day  24  is  consistent  with  the 
poo.oiiity  that  gradual  accumulation  of  matrix  down-regulates  matrix  biosynthesis. 
While  biosynthetic  rates  were  similar  for  chondrocytes  in  agarose  and  cartilage,  the 
GAG  accumulation  rate  per  disk  volume  in  agarose  (Table  1),  calculated  by  linear 
regression  of  GAG  content  over  the  10-week  period  (Fig.  5),  was  lower  than  that  in 
intact  cartilage.  The  rate  at  which  GAG  was  released  from  the  agarose  disks  was  also 
lower  1  Table  1  i.  These  results  are  consistent  with  the  fact  that  the  cell  density  in  the 
agarose  disks  is  still  4  5  times  iower  than  that  in  the  parent  cartilage  tissue. 


1 1 


GAG  (yg/plug) 
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A  CONTENT 


GAG/ 


DNA 


o  i  - 1 - 1 - 1 - 1 - 1 - 1 - 1  o.5 

10  20  30  40  50  60  70 

Day  of  Culture 


Figure  3.  Chondrocytes  in  agarose  :  GAG  and  DNA  content  during  culture. 
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Figure  4.  Chondrocytes  in  agarose  :  geometrical  properties  during  culture. 
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Thickness 


Biosynthesis  Rate 


Day  of  Culture 


Figure  5.  Chondrocytes  in  agarose  :  radiolabel  incorporation  during  culture. 


Cell-laden  disks  showed  significantly  enhanced  equilibrium  modulus,  dynamic 
stiffness,  and  streaming  potential  with  increasing  time  in  culture  compared  to  con¬ 
trols  (Figs.  6  and  7),  consistent  with  the  increased  GAG  content  found  in  the  same 
specimens.  The  existence  of  the  streaming  potential  and  the  increase  in  the  poten¬ 
tial  with  time  in  culture  is  especially  significant,  since  this  streaming  potential  gives 
a  nondestructive  measure  of  the  quality  of  the  newly  synthesized  extracellular  ma¬ 
trix  secreted  and  assembled  by  the  cells.  If  the  newly  synthesized  proteoglycans  had 
not  been  functionally  immobilized  within  the  developing  matrix  (as  occurs  in  intact 
cartilage),  there  would  have  been  little  or  no  increase  in  streaming  potential  (Fig.  7). 

The  frequency  response  of  the  streaming  potential  and  dynamic  'riffness  is 
also  an  important  measure  of  the  functional  behavior  of  the  extracellular  matrix.  The 
dynamic  stiffness  and  streaming  potential  of  agarose/chondrocyte  disks  increased  with 
frequency  (Fig.  8)  in  a  manner  characteristic  of  the  poroelastic  behavior  of  cartilage 
37,  23] .  Sinusoidal  stress-strain  and  streaming  potential  were  essentially  linear  in  the 
0.01-1  Hz  range  (total  harmonic  distortion  <  8%).  The  frequency  dependence  of  the 
dynamic  stiffness  and  streaming  potential  amplitude  of  a  disk  of  normal  cartilage  is 
also  shown  in  Fig.  8  for  comparison.  The  increase  in  streaming  potential  amplitude 
with  increasing  frequency  is  due  to  the  increase  in  relative  fluid  velocity  produced  at 
higher  compression  rates.  This  relative  fluid  velocity  is  responsible  for  producing  the 
'teaming  potential. 

The  equilibrium  modulus  and  the  computed  hydraulic  permeability  of 
agarose/ chondrocyte  disks  at  day  50  are  compared  to  values  for  parent  calf  cartilage 
disks  in  Table  1.  The  permeability  is  slightly  higher,  consistent  with  the  larger  pore 
size  and  lower  GAG  content  of  the  chondrocyte/ agarose  disk  compared  to  cartilage. 
Swarm  Rat  Chondrosarcoma  Cells 

In  another  set  of  experiments,  Swarm  rat  chondrosarcoma  cells  incor¬ 

porated  into  2%  agarose  gels  as  before.  During  4  weeks  in  culture,  the  total  GAG 
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Day  of  Culture 


Figure  6.  Chondrocytes  in  agarose  :  stiffness  during  culture. 
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Figure  7.  Chondrocytes  in  agarose  :  streaming  potential  during  culture. 
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POTENTIAL,  nV/%  STIFFNESS  ,  MPa 


Figure  8.  Chondrocytes  in  agarose  :  stiffness  and  streaming  potential  of  chondro¬ 
cyte  agarose  plugs  displays  poroelastic  behavior  like  articular  cartilage. 
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Table  2.  CHONDROCYTES  IN  AGAROSE  :  PHYSICAL  PROPERTIES  ON  DAY  50 


COMPARISON  : 


CHONDROCYTE/AGAROSE  DISKS 
and  CALF  CARTILAGE  EXPLANTS 


Chondrocyte/  Calf  Cartilage 
Agarose  diskst  'parent  disks’ 

!  Equilibrium  Modulus  (kPa)  91  ±  7  (3)  470  [1 ,2] 

Dynamic  Stiffness  (0.5Hz)  (MPa)  1.64  ±  0.31  (3)  4.7  [1,2] 

^Streaming  Potential  (0.5Hz)  (MV/%)  61  ±  17(3)  150-600  [1,2] 

Hydraulic  Perm  (m4/N  sx  1015)  5.1  ±1.9  (3)  3.5  [1,2] 

Electrokinetic  Coefficient  (mV/MPa)  3.7  ±  .8  (3)  20  [1 ,2] 


t  data  for  day  50 


MEAN  ±  SD  (n) 


1)  Frank  EH  &:  Grodzinsky  AJ.  J  Biomech,  20:629,  1987. 

2)  Hey  LA  BS  Thesis,  MIT,  1984. 


content  within  the  agarose/chondrocyte  disks  increased  almost  20  fold,  while  DNA 
content  increased  by  2-3  fold  (data  not  shown).  After  4  weeks,  dynamic  stiffness  and 
streaming  potential  of  cell-laden  disks  had  also  increased  above  control  levels.  These 
experiments  were  actually  performed  before  those  involving  normal  calf  chondrocytes. 
Because  of  changes  and  improvements  in  the  methodology  in  mixing  and  casting  the 
gels,  it  is  difficult  to  compare  the  absolute  value  of  changes  observed  with  normal 
chondrocytes  and  chondrosarcoma  cells  thus  far.  It  is  clear,  however,  that  definitive 
changes  have  been  observed  with  both  types  of  cells. 

Conclusions 

The  results  of  these  studies  suggest  that:  (1)  both  normal  chondrocytes  and 
Swarm  rat  chondrosarcoma  cells  in  agarose  culture  can,  under  proper  culture  con¬ 
ditions,  continue  to  synthesize  matrix  macromolecules  at  a  rate  similar  to  that  in 
native  cartilage,  and  (2)  chondrocytes  in  agarose  can  successfully  mediate  the  as¬ 
sembly  and  accumulation  of  a  normal,  mechanically  functional  extracellular  matrix. 
These  results  provide  support  for  the  long-term  objectives  of  this  research  in  that 


the  chondrocyte/agarose  system  appears  to  be  an  appropriate  model  with  which  to 
study  the  mechanisms  by  which  physical  forces  (e.g.,  electrical  and  mechanical)  may 
interact  with  mammalian  cells.  The  chondrocytes  can  be  maintained  in  a  natural 
environment  which,  in  certain  ways,  is  more  amenable  than  that  of  intact  cartilage 
to  the  study  of  mechanism  at  the  single  cell  level. 

MODEL  OF  EQUILIBRIUM  PROPERTIES 

Summary 

A  microcontinuum  model  based  on  the  unit-cell  model  of  rod-like  polyelec¬ 
trolytes,  and  representing  a  chondroitin  sulfate  molecule,  is  used  to  predict  osmotic 
swelling  pressures  in  proteoglycan  solutions  and  the  ionic  strength  dependent  compo¬ 
nent  of  the  equilibrium  modulus  of  articular  cartilage.  The  model  predictions  are  com¬ 
pared  to  Donnan  theory  for  proteoglycan  solutions.  Donnan  theory  predicts  higher 
pressures  and  is  inconsistent  with  the  rod-model  prediction.  Donnan  theory  may  be 
made  to  fit  experimental  data  by  ad-hoc  usage  of  activity  coefficients  and  unphysical 
designation  of  excluded  volume  pressure.  The  rod-model  only  requires  knowledge  of 
nonthermodynamic,  structural  parameters  which  may  be  measured  independently  of 
the  theory.  The  rod-model  accurately  describes  the  osmotic  pressure  of  proteoglycan 
solutions  and  the  ionic  strength  dependent  component  of  the  equilibrium  modulus 
of  articular  cartilage  with  parameters  which  are  reasonable  when  compared  to  those 
of  chemical-structural  models  of  chondroitin  sulfate,  or  those  obtained  from  other 
experiments. 


Introduction 

The  model  considered  here  is  the  mirrornn<inuum  unit  cell  [39,  40,  41]  where 
the  solid  is  represented  as  a  circular  cylinder  with  a  uniform  surface  charge  and  the 
electrolyte  fluid  surrounds  the  solid  out  to  a  fictitious  cylindrical  boundary,  the  radius 
of  which  is  determined  by  the  solid  to  fluid  volume  ratio  (Fig.  9).  The  osmotic  swelling 
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Figure  9.  Unit  Cell. 

pressure  and  ion  distributions  will  be  calculated  using  the  Poisson- Boltzmann  (PB) 
equation.  The  PB  equation  is  an  approximation  to  the  exact  many  body  problem 
42  .  In  the  PB  approximation  the  solvent  is  treated  as  an  incompressible  dielectric 
continuum  and  the  potential  of  mean  force  on  an  ion  is  equated  with  the  electro¬ 
static  potential.  The  PB  equation  excludes  effects  such  as  dielectric  saturation  and 
electrostriction  while  the  more  exact  model  includes  ionic  polarizability,  short-range 
interionic  repulsion,  and  any  other  2-particle  correlations  beyond  the  mean  field, 
point  ion  approximation.  However,  many  theoretical  and  experimental  studies  have 
confirmed  the  accuracy  of  the  PB  equation  in  a  variety  of  situations  [43]. 

Comparison  to  Monte  Carlo  simulations  where  the  electrolyte  model  the  sim¬ 
ple  primitive  model  (incorporates  hard  sphere  ion-ion  repulsion)  have  suggested  that 
the  PB  approximation  is  valid  when  the  concentration  of  divalent  ions  is  low  [41,  44]. 
Divalent  ions  have  more  significant  2-particle  correlations  which  push  them  closer  to 
the  charged  surface  than  is  predicted  by  the  PB  equation,  reducing  the  osmotic  pres¬ 
sure.  Experiments  measuring  surface  potentials  and  forces  between  charged  surfaces 
have  agreed  well  with  PB  theory  [43]. 
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Condensation  theory,  another  approach  to  the  properties  of  polyelectrolytes, 
has  been  formulated  using  empirical  rules  by  Manning  [30]  and  Oosawa  [45].  It  has 
been  shown  to  be  a  special  case  of  the  PB  formulation  and  n“ed  not  be  used  since  the 
PB  equation  is  not  difficult  to  solve  numerically  and  the  PB  approximatioii  is  more 
general  in  consequence  than  Condensation  theory  [46,  47]. 

Numerical  Solution  of  the  Poisson-Boltzmann  Equation 

The  electrostatic  potential  <j>  and  charge  density  p  satisfy  Poisson’s  equation 

~  (i) 

with  the  medium  permittivity  e  .The  *  represents  a  dimensional  (unnormalized) 
quantity.  The  charge  in  the  fluid  is  mainly  due  to  dissociated  sodium  chloride 
(NaCl)  and  is  thus  represented  as  satisfying  the  Boltzmann  distribution  for  a  mono- 
monovalent  electrolyte. 

p  =  F(C+  -  C_) 

=  FCo(eT&  -  e*&)  (2) 

=  -2FC0  sinh(^|) 

where  F  is  the  faraday  constant,  R  the  gas  constant,  T  temperature,  and  Co  the 
concentration  of  NaCl  in  the  bath  where  <j>  is  set  to  zero.  The  PB  equation  is  then 
obtained. 

H££2sinh(|i)  (3) 

The  boundary  condition  on  the  solid  cylinder  surface  ( r  =  a)  arising  from  Gauss’s 
law  is 

&  =  eo  +•£*.«♦  *  (m-Fr,a-  (4) 

where  a  is  the  uniform  surface  charge,  Er  is  the  r&diai  component  of  the  electric  field, 
a+  indicates  just  outside  the  cylinder,  and  a-  just  inside.  Inside  the  cylinder  the  field 


22 


is  laplacian,  and  is  given  by 


4>in  =  A  +  ]T](f?nfn  cos  nO  +  Cnfn  sin n.0).  (5) 

n  =  l 

Due  to  the  azimuthal  symmetry,  all  B’s  and  C’s  are  zero  so  that  there  is  no  field 
inside.  Then, 


a 


Za+  E rta+ 

d<}>\ 


—  £ 


dr 


At  the  outer  cell  boundary  (f  =  &), 


dij> 

dr  g 


=  0 


(6) 

(7) 


(a) 


is  required  by  considering  two  unit  cells  butted  together  (without  equation  (8)  there 
would  be  a  sheet  of  charge  surrounding  each  unit  cell).  With  the  normalization 


where 


<7  = 
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RTe 
2  C0F* 


<?o 


eRT 

~w 


and  where  8  is  the  debye  length,  equations  (3),  (7),  and  (8)  become 


(9) 

(10) 

(11) 


(12) 

(13) 


<P<f>  1  dd 

dr 3  r  dr 


sinh(d>) 


(14) 
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—  a 


(15) 


dd> 

Tr^ 0  (16) 

There  is  no  known  general  analytic  solution  to  equation  (14)  [48]  and  numerical 
solutions  (Runge-Kutta)  can  be  nonconvergent  in  the  nonlinear  regime,  even  with 
very  small  steps  [49]. 

An  efficient  numerical  solution  can  be  formulated  by  solving  equation  (14)  in 
intervals  where  the  change  in  <p  <C  1.  Considering  an  interval  (ro,^),  equation  (14) 


is  linearized  to 


where  <f>0  =  4>{r o)-  Then, 


=  <t>  -  <Po 


sinh($)  =  sinh(-0)  cosh(<^o)  +  cosh(V>)  sinh(^o) 


=  ip  cosh(^o)  +  sinh(<^o) 


and  equation  (14)  becomes 


with  the  solution 


and  derivative 


d?ip  1  dip 

tt  H - j - ipcosh(ipo)  =  sinh((/»o) 

ar*  t  dr 


ip  =  Cilo(ar)  +  C3K0(ar)  -  tanh(<^0) 


=  aCiIi(ar)  —  aC2Ki(ar) 


where  In  and  Kn  are  modified  Bessel  functions  of  order  n,  and 


a2  =  cosh(0o)- 
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C\  and  C2  can  be  found  with  knowledge  of  <j>0  and  4>'0.  So  that  at  r0, 


V>(ro)  =  o 

(24) 

di’(ro)  ±, 

dr  =*<> 

(25) 

become 

Cil0(ar0 )  4-  C2K0(arQ)  =  tanh(^0) 

(26) 

Clh{ar0)-C2K1(ar0)  =  ^ 

a 

(27) 

which,  solved  for  Cj 

and  C2  is 

<t>'a  ,  K\  (orro  )  t&nh(^o  ) 

a  '  Ko(aro) 

'  “  «<*■„)  + 

(28) 

tanh(<^o)  ~  Cil0(ar0) 

2  ~  Ko(ar0) 

(29) 

A  step  is  taken  to  the  new  point  by  calculating  ^(rx)  and  from  equations 

(21)  and  (22),  which  then  serves  as  the  starting  point  r0  for  the  next  step,  with  a 
restriction  on  the  nonlinearity  (second  order  term  in  (19)), 


2 


< 


(30) 


If  the  trial  step  violates  equation  (30),  the  step  size  is  reduced  until  equation  (30) 
is  satisfied.  The  complete  integration  begins  at  the  outer  boundary  (r=b)  where 
(16)  specifies  <f>'0  —  0  and  a  guess  is  taken  for  <f>0  —  <j>{b).  Equation  (15)  is  satisfied 
by  repeating  the  integration  with  different  values  of  this  initial  <f>(b)  determined  by 
polynomial  interpolation  of  the  pairs  (<^>(5),  )  until  =  —a  to  within  a  specified 

accuracy  (~  10-6).  By  using  a  five  point  discrete  differentiation  formula  to  calculate 
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the  error  in  the  solution  to  equation  (14), 


error  — 


+  lg-sinh(j>) 
sinh(<£) 


(31) 


the  error  was  found  to  be  of  the  same  order  as  which  was  chosen  to  be  10  6. 


Electrostatic  Swelling  Pressure  in  the  Unit  Cell 

It  is  proposed  that  the  swelling  pressure  due  to  the  electrostatic  repulsion  of 
an  array  of  charged  rods  is  the  pressure  (osmotic)  at  the  cell  boundary  minus  that  in 
the  bath, 


P.Weii  =  RT{c+(b)  +  c_(fe)  —  2C0)  (32) 


=  2RTC0(cosh{<f>(b)  -  1)) 


(33) 


This  result  has  been  shown  by  differentiating  the  free  energy  (comprised  of  the  elec¬ 
trostatic  energy  of  all  the  charge  in  the  unit  cell  and  the  entropy  of  mixing  of  the 
mobile  charge)  with  respect  to  a  change  in  volume  of  the  unit  cell  [40].  Here,  the 
result  of  equation  (33)  will  be  found  from  the  consideration  of  the  fluid  stress  tensor 
in  the  unit  cell  with  no  reference  to  a  cell  free  energy  or  hypothetical  charging  process. 
The  two  perspectives  are  complementary. 


jj-i  =  +  ViP  +  RT\n(~fiXi)  +  ZiF4>  (34) 


where  =  ^5^-  (rij  =  #  moles  of  species  i)  is  the  mole  fraction,  Vi  =  ^~\T p  ^ 
the  partial  molar  volume,  P  the  pressure  of  the  mixture,  7<  the  activity  coefficient,  z,- 
the  valence,  and  <j>  the  electrostatic  potential.  Neglecting  temperature  gradients  and 
assuming  ideality  (7^  =  1),  the  chemical  potential  of  the  uncharged  solvent  is, 

£0  =  v0P  +  RTln  x0.  (35) 
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By  definition, 

1  =  x0  +  X+  +  z_, 

(36) 

and  for  dilute  solutions, 

X+  +  X-  <  1, 

(37) 

so  that  one  gets 

In  x0  =  — (x+  +  x_). 

(38) 

Then, 

p.0  =  v0P  -  RT(x+  +  z_) 

(39) 

For  small  solutes  (NaCl)  the 

ViP  term  may  be  neglected  [50], 

=  RTlnx++F<i> 

(40) 

ji _  =  RT  In  z_  -  F4>. 

(41) 

Each  species  is  in  equilibrium  in  the  unit  cell  and  between  the  cell  and  the  bath 

(through  the  array  of  other  rods)  requiring, 

VA,  =  0 

(42) 

Ai  =  constant  =  0 

(43) 

since  the  constant  is  arbitrary.  Then  equation  (39)  defines  the  pressure, 


p  =  RT(X+  +  J-) 

v0 


RT(c+  +  c_) 


(44) 


The  equilibrium  force  balance  may  be  found  from  equations  (40),  (41),  and  (42), 


0  =  —  c+  V/z_  —  c'+  V/i  _ 

=  -12TV(x+  +  r.)^(r+-i.)(-% 


27 


(45) 

(46) 


Dividing  by  v0, 

0  =  -RTV(c+  +  c'_)  +  F(c+  -  cl)(-V<£).  (47) 

Using  equation  (44)  along  with, 

F{c\  -  c_)(-V<£)  =  pE  =  V-7  (48) 

where  T  is  the  electrostatic  stress  tensor, 

7  =  eEE  -  ^lE2  (49) 

and  I  the  idemfactor,  equation  (47)  becomes 

0  =  V  •  {  —  Pi  +  7),  (50) 

which  is  the  first  integral  of  equation  (3),  the  PB  equation  (multiply  equation  (3)  by 
Vd>).  The  total  stress  tensor  is  then, 

F  =-Pl  +  f.  (51) 

The  force  on  any  volume  may  be  found  by, 

F  =  fsT-dS  (52) 

where  5  is  the  enclosing  surface.  There  is  no  net  force  on  the  solid  rod  due  to 
the  azimuthal  symmetry.  This  can  be  seen  by  constructing  a  volume  enclosing  the 
rod  at  its  surface.  In  an  actual  array  of  cylindrical  rods,  only  the  central  rod  will 
possess  this  symmetry  and  experience  no  net  force.  For  non-central  rods  there  will 
be  an  asymmetry  which  will  give  rise  to  a  net  force  on  the  rod  (pushing  it  outwards) 
requiring  a  mechanical  restraint  (like  the  solid  structure). 
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By  forming  a  surface  on  the  outer  boundary  of  the  unit  cell  where  7  =  0  (due 
to  equation  (8))  it  is  seen  that  although  no  part  of  the  unit  cell  experiences  any  net 
force,  there  is  a  pressure 


h  =  +  c_(S)) 

-  2RTC0  cosh(<£(6)) 

which  is  higher  than  the  pressure  in  the  bath, 


P0  =  2  RTC0. 


(53) 

(54) 


(55) 


And  it  is  the  difference, 

P.W'U  =  A  -  Po  =  2Ai  c-'o(cosh(0(6)  -  1))  (56) 

which  constitutes  the  swelling  pressure  that  is  transmitted  to  the  solid  array  by  the 
previously  mentioned  asymmetries  . 

Comparison  to  Chondroitin  Sulfate  Solutions 

The  swelling  pressure  of  chondroitin  sulfate  solutions  has  been  measured  by 
Urban  et  al.  [27].  Commercially  obtained  CS  preparations  (purity  ~  70%)  were 
placed  in  a  dialysis  sac  and  equilibrated  with  a  Polyethylene  Glycol  (PEG)  solution 
of  known  osmotic  pressure.  Radioactive  tracer  ions  were  used  to  measure  the  fixed 
charge  density  (  based  on  the  Donnan  model  ).  The  swelling  pressure  of  the  solution 
as  a  function  of  fixed  charge  density  was  then  compared  with  that  predicted  by  a 
macroscopic  ideal  Donnan  equilibrium  model  at  2  ionic  strengths,  0.15  M  and  1.5  M. 
The  expression  for  swelling  pressure  in  the  ideal  Donnan  model  follows  from  bulk 
electroneutrality  and  ideal  Donnan  equilibrium, 

Pm  +  Cf  -  CD  =  0  (57) 
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(58) 


=  cl 

where  pm  is  the  macroscopic  fixed  charge  in  moles  per  litre  tissue  fluid,  C±  are  the 
Donnan  ion  concentrations  in  the  polyelectrolyte  phase  and  C0  is  the  concentration 
in  the  equilibrating  bath.  C±  are  found  from  equations  (57)  and  (58)  , 

Cl  =  W2  +  y/fZJA  +  cj  (59) 

from  which  the  swelling  pressure  is  taken  as  the  osmotic  pressure  in  the  polyelectrolyte 
phase  minus  that  in  the  bath, 


Pswell  =  RT{Cl  +  CD_  -  2 Ca) 

(60) 

=  =  RT(2yJf,in  +  c;  -  cl) 

(61) 

in  units  of  kPa  (  101.325  kPa/atm  ).  The  data  is  shown  along  with  the  ideal  Donnan 
prediction  in  Figures  10  and  11. 

According  to  the  macroscopic  model,  the  discrepancy  between  model  predic¬ 
tion  and  data  in  Figure  10  can  be  accounted  for  by  introducing  activity  coefficients 
in  an  ad-hoc  manner.  Wells’  modifications  of  Manning’s  expressions  [27,  28)  predict 
a  swelling  pressure  much  lower  than  that  measured. 

The  underprediction  of  the  ideal  Donnan  model  in  Figure  11  is  suggested  [27] 
to  be  due  to  excluded  volume  effects.  However,  at  1.5  M  with  a  debye  length  of  ~  1.4 
A,  virtually  no  mean  field,  point  ion,  electrostatic  contributions  to  swelling  pressure 
will  exist.  Then  the  total  swelling  pressure  in  Figure  11  is  interpreted  to  be  due  to 
excluded  volume,  and  other  forces  not  accounted  for  in  the  PB  approximation;  it  is 
the  difference  in  swelling  pressure  at  0.15  M  and  1.5  M  that  is  due  to  electrostatic 
repulsion  and  predicted  by  the  unit  cell  model. 

Figure  12  shows  both  data  sets  with  polynomial  best  fits.  The  subtraction 
of  the  2  curves  at  4  points  is  plotted  in  Figure  13  along  with  a  best  fit  micromodel 
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Pressure,  atm 


Swelling  Pressure  of  CS  vs.FCD 
at  0.15  M  NaCI 


•  Fixed  Charge  Density,  meq/gmH20 

Urban  JPQ  Maroudaa  A  Blorhaology  vol.  IS  pg  4S7  fig  5, 1979 


Figure  10.  Swelling  pressure  of  chondroitin  sulfate  solution  in  0.15  M  NaCI. 


31 


Pressure,  atm 


Swelling  Pressure  of  CS  vs.  FCD 
at  1.5  M  NaCI 


Urban  JM  Maroodaa  A  Blortwotogy  vol.  16  pg  497  fig  6, 1976 


Figure  11.  Swelling  pressure  of  chondroitin  sulfate  solution  in  1.5  M  NaCI. 
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Table  3.  MICROMODEL  FIT  (on  a)  WITH  a  =  6  A  FIXED 


a  = 

6  A 

_ 

-0.1379  C/m2 

chi  =  0.04569 

ma 

b 

0 b-a ) 

model  P 

data  P 

0.35 

22.92 

16.92 

0.6017 

0.5744 

0.45 

20.41 

14.41 

1.283 

1.159 

0.55 

18.64 

12.64 

2.189 

2.086 

0.65 

17.31 

11.31 

3.276 

3.414 

prediction.  The  parameters  required  to  solve  equation  (3)  to  obtain  the  swelling 
pressure  are  the  radii  a,  b,  the  surface  charge  ir,  and  the  ionic  strength  C0.  The 
permittivity  c  is  that  of  water  at  25 °C , 

e  =  78.Zea  =  6.93  x  10~10  Farad/m  (62) 

The  unspecified  parameters  in  the  unit  cell  model  are  then  a,  b,  and  a,  subject  to  the 
constraint, 

.  _  2  era 

Pm  ~  103F(62  -  a2)’ 

leaving  2  degrees  of  freedom.  The  variable  parameters  were  chosen  to  to  be  a  and 
1 7  so  that  b  is  determined  from  equation  (63).  The  model  was  then  fit  using  a  2 
dimensional  minimization  of  the  sum  of  the  squares  of  the  differences  (chi)  between 
data  and  model  at  4  points  on  the  —  pm  axis.  The  data  could  be  fit  with  more  than 
one  set  of  (a,  a),  depending  on  the  choice  of  starting  point  and  end  condition.  By 
fixing  the  rod  radius  a  and  performing  a  1-dimensional  minimization  (on  a),  the  3 
best  fits  in  Tables  3,  4,  5  were  found.  It  is  the  fit  in  Table  3  that  is  shown  in  Figure 
13. 

As  there  are  no  standard  deviations  on  the  data,  confidence  limits  for  the  fits 
cannot  be  stated.  It  is  most  likely  that  none  of  the  3  shown  actually  better  represents 
the  data  than  the  other  two.  Thus,  apparently  there  is  really  only  1  degree  of 


Pressure,  ilm 


Swelling  Pressure  of  CS  vs  FCD 


•  Fixed  Charge  Density,  meq/gmH20 


Figure  12.  Polynomial  fits  to  swelling  pressure. 
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Pressure,  atm 


Swelling  Pressure  of  CS  vs  FCD 
Micromodel  Fit 


Rod  Radius  >6A 


Figure  13.  Micromodel  fit  to  the  electrostatic  part  of  the  swelling  pressure. 
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freedom  in  the  fit  because  variations  in  the  2d  degree  of  freedom  would  not  produce 
a  significantly  better  chi  (if  standard  deviations  were  available).  Had  other  data 
been  reported  (such  as  hydration,  which  the  procedures  indicate  was  taken)  another 
parameter  (ex.  d)  could  have  been  fixed  by  the  experiment. 

Here  since  there  is  effectively  only  1  degree  of  freedom  in  the  fit,  the  other  must 
be  chosen.  Literature  values  indicate  a  hydrodynamic  radius  of  5  — 8.5  A°  [51],  whereas 
)  a  value  9  A°  was  found  to  agree  with  measurements  of  hydraulic  permeability  and 

electrokinetic  coefficients  of  cartilage  [24].  Tables  3,4,  and  5  show  that  the  data  may 
be  fit  (by  varying  a)  with  a  satisfying  any  of  these  suggested  values.  It  is  pertinent 
to  note  that  the  quantity  (6  —  a)  was  nearly  constant  for  the  different  fits  suggesting 
that,  as  far  as  swelling  pressure  is  concerned,  it  is  the  separation  distance  of  the 
charged  surfaces  that  is  most  critical  (not  a  or  b  separately),  which  could  explain  the 
degeneration  from  2  to  1  degree  of  freedom. 

Comparison  to  the  Equilibrium  Modulus  of  Articular  Cartilage 

The  ionic  strength  dependence  of  the  equilibrium  modulus  of  articular  car¬ 
tilage  is  a  manifestation  of  long-range  electrostatic  interactions  [22].  At  high  ionic 
strength  (>  1 M)  where  long-range  electrostatic  interactions  are  shielded,  the  equilib¬ 
rium  modulus  can  be  considered  as  arising  from  forces  not  accounted  for  in  the  unit 
ceil  PB  approximation.  As  the  ionic  strength  is  lowered,  the  increase  in  modulus  is 
interpreted  as  a  result  of  the  repulsion  of  fixed  charge  groups  on  the  matrix.  At  phys¬ 
iological  ionic  strength  (0.15A/)  the  modulus  is  about  twice  that  found  at  1 M  (0.55 
MPa,  0.27  MPa).  With  the  assumption  that  the  microstructure  of  the  matrix  does  not 
change  significantly  as  ionic  strength  is  changed  and  can  be  adequately  represented 
as  the  same  unit  cell,  the  ionic  strength  dependent  contribution  to  the  equilibrium 
modulus  can  be  modelled  as  due  to  change*  in  the  swelling  pressure  predicted  by 
equation  (33)  in  response  to  strain.  Then 

H\  =  Ha  -  H ?  (64) 
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where  H A  is  the  equilibrium  modulus,  H\  the  electrostatic  component,  and  H *  the 
equilibrium  modulus  at  high  ionic  strength.  In  this  section  the  model  will  be  fit  to 
the  data  of  Figure  8,  pg  156  of  Eisenberg  and  Grodzinsky  [22].  This  data  itself  is  a 
fit  to  descriptive  functions  not  representing  any  model.  The  standard  deviations  are 
the  deviations  from  this  fit  for  1  to  5  cartilage  plugs  6.4  mm  in  diameter  taken  from 
the  femoropatellar  region  of  2-year-old  cattle. 

The  testing  was  done  in  a  range  of  strain,  0%  to  25%,  so  that  the  modulus 
from  the  unit  cell  is  found  here  as 


P,well(-2)  ~  1) 

0.1 


(65) 


where  P,weu(. 2)  is  the  swelling  pressure  at  20%  strain,  P,weu(.  1)  at  10%  strain.  One 
could  choose  different  and  multiple  strains  to  find  HeA,  but  the  result  would  differ 
little.  Equation! 65 )gives  a  modulus  (independent  of  strain)  valid  in  the  strain  range 
experimentally  tested. 

To  find  H\,  one  must  choose  a,  b„  (the  unstrained  cell  boundary),  and  a  as 
well  as  a  model  for  strain,  a  is  in  the  range  5  — ►  8.5/1°  [51].  Given  the  solid  volume 
fraction  of  the  GAG  constituent  V,  (Appendix)  b0  can  be  found, 


<7  will  be  the  fitted  parameter.  The  strain  will  be  modeled  as  shortening  the  outer  cell 
diameter  to  b  from  the  unstrained  value  b0.  Given  a  strain  eG  on  the  GAG  constituent, 


&  =  60(1-sg)1/2. 


(67) 


The  GAG  strain  is  related  to  total  tissue  strain  t  by 


eG  =  ke 


(68) 
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where 


1  -  Vc 


and  \VT  is  the  fractional  water  content  of  the  tissue  as  a  whole  (Appendix). 

Choosing  a  rod  radius  of  6  A,  a  solid  volume  fraction  of  0.09,  and  Wj  —  0.8 
the  best  fit  in  Figure  14  was  found  by  minimizing  chi  with  respect  to  variations  in  cr. 
The  surface  charge  corresponds  well  with  structural  models  of  the  chondroitin  sulfate 
chain.  The  structural  model  predicts  an  intercharge  spacing  of  0.5 nm  —*  0.7am 
51,  52!  whereas  the  surface  charge  found  by  the  fit  corresponds  to  spacing  of  0.73nm. 
The  reasonable  fit  between  theory  and  experiment  suggests  that  electrostatic  repulsion 
at  the  GAG  chain  level  (as  predicted  by  the  PB  equation)  is  the  basis  of  the  ionic 
strength  dependent  component  of  the  equilibrium  modulus. 


EFFECT  OF  APPLIED  ELECTRIC  FIELDS 


Experimental 

Given  the  demonstrated  development  of  a  mechanically  functional  extracel¬ 
lular  matrix,  we  have  now  begun  a  series  of  experiments  to  study  the  transduction 
mechanisms  involved  in  response  to  electrical  stimulation. 

Initial  studies  have  utilizedTeflon  chambers  (Fig.  15  )  each  of  which  contains 
two  independent  lanes  capable  of  holding  16  cylindrical  agarose/chondrocyte  disks  in 
each  lane.  Gel  disks  are  placed  intoTeflon  disk  holders  which,  in  turn,  are  inserted  into 
the  lanes  as  shown  in  Fig.  15.  Two  prototype  chambers  have  now  been  constructed, 
but  more  will  be  needed.  The  disk  holders  will  accept  3mm  diameter  by  1mm  thick 
disks.  (These  chambers  were  constructed  as  part  of  a  Bachelor’s  Thesis  project  by 
Paul  M.  Tiao,  June  1989.)  The  lanes  in  the  chambers  are  connected  to  platinum  elec¬ 
trodes  by  means  of  “salt  bridges”  comprised  of  autoclavible  plastic  tubing  filled  with 
5%  w/v  agarose.  The  electrode  baths  are  filled  with  saline.  This  external  electrode 
arrangement  isolates  the  agarose/chondrocyte  disks  from  any  electrode  reaction  prod¬ 
ucts.  (Our  previous  experiments  using  cartilage  organ  culture  explants  showed  that 
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Modulus  (ES  comp)  vs.  Ionic  Strength 
Micromodel  Fit 


Rod  Radius  =  6  A 
Call  Radius  *  20  A 
Surfac#  Chargt  =  <0.0580  C/m2 
Chi  s  2.51 


Figure  14.  Electrostatic  component  of  articular  cartilage  equilibrium  modulus  micro¬ 
model  fit  to  data  of  Figure  8,  pg  156  of  Eisenberg  and  Grodzinsky  1985. 
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Figure  15.  Electrical  Stimulation  Chamber. 


such  electrode  reaction  products  could  themselves  stimulate  the  synthesis  of  stress 
response  proteins  in  the  cartilage.) 

Further  work  will  have  to  be  done  to  quantify  temperature  at  various  loca¬ 
tions  within  the  stimulation  chamber.  Previous  tests  of  this  system  have  shown  that 
negligible  ohmic  heating  occurs  for  current  densities  up  to  10-20  mA/cm2;  feedback- 
controlled  circulating  cooling  fluid  may  be  needed  for  higher  current  densities. 

Since  it  is  absolutely  essential  that  temperature,  field  and  current  parameters 
be  continuously  monitored  and  recorded  throughout  each  test,  a  dedicated  microcom¬ 
puter  system  will  have  to  be  developed  for  this  purpose. 

A  preliminary  experiment  on  the  effect  of  electrical  currents  on  sulfate  and 
proline  incorporation  for  day  32  plugs  is  shown  in  Fig.  15. 

While  radiolabel  incorporation  appears  to  change  with  applied  current  density, 
additional  experiments  will  have  to  be  performed  to  confirm  this  result. 
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Figure  16.  Electrical  Stimulation  at  100  Hz  of  3mm  Chondrocyte/ Agarose  plugs  from 
day  32  and  its  Effect  on  Proline  and  Sulfate  Incorporation. 
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APPENDIX 


CARTILAGE  PARAMETERS 

Cartilage  Constituent  Compartments 

To  implement  the  unit  cell  to  describe  cartilage  properties,  an  account  must 
be  made  of  the  various  constituents  of  cartilage,  namely  the  GAG  constituent  as 
represented  bv  the  unit  cell,  and  everything  else  (collagen  etc.).  There  is  both  solid 
and  fluid  in  each  giving  rise  to  four  distinguished  volume  fractions  : 


SG  - 

gag  solid 

Wo  - 

gag  water 

ho  - 

other  solid 

\V0 

other  water 

Only  3  equations  are  independent  since  they  sum  to  1.  Combinations  of  the  3  are 
available  from  the  literature, 


WT  =  Wa  +  Wo  =  .70  .80 


(70) 


53.  54 


,  \V0 

JOW  "  Tvv  =  35 


(71) 


fcs  - 


’<7 


So  d  Sc, 


.15  -a  .25. 


(72) 


For  model  implementation  the  following  3  parameters  (and  their  ranges  from  above) 
are  most  useful, 

WT  =  .70  .80  (73) 

WG 


few  =  ppm  =  1  -  fow  =  .6  .7 


(74) 
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(75) 


The  strain  in  the  GAG  constituent  is  more  than  that  of  the  total  tissue  due  to 
its  higher  water  content.  With  the  assumption  that  / gw  (or  fow)  is  constant  while 
fluid  is  lost  from  the  tissue  during  strain, 

£g  _  1  ~  Vs 
e  ~  WT 

Where  e  is  the  total  tissue  strain  and  sg  is  the  strain  on  the  GAG  constituent.  For 
example,  when  all  the  fluid  is  excluded,  e  =  Wj  ~  0.8  and  cq  =  1  —  Vs  ~  0.92  and 
only  the  solid  remains. 


52 


*  U.  S.  GOVERNMENT  PRINTING  OFFICE:  1 99 1 --56 1 -052/4 0004 


